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Abstract 

We derive, via simultaneous homogenization and dimension reduction, 
the F-limit for thin elastic plates whose energy density oscillates on a scale 
that is either comparable to, or much smaller than, the film thickness. 
We consider the energy scaling that corresponds to Kirchhoff 's nonlinear 
bending theory of plates. 
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1 Introduction 

KirchhofF's nonlinear plate theory associates with a deformation u : 5 — > K'^ of 
a two-dimensional stress-free reference configuration 5 C the bending energy 

Q2mx'))dx', (1) 

s 

where Q2 is the quadratic form from linear elasticity, and II denotes the sec- 
ond fundamental form associated with u. The key condition on the admissible 
deformations u is that they must satisfy the isometry constraint 

dau- dpu = 5ai3, a,/3e{l,2} (2) 

where Sap denotes the Kronecker delta. Physically, (1) describes the elastic 
energy stored in a deformed plate that can undergo large deformations but not 
shearing or stretching. 

In [FJM02] Kirchhoff 's nonlinear plate theory was rigorously derived as a zero- 
thickness T-limit from 3d nonlinear elasticity. In this article we combine their 
result with homogenization. We consider a plate with thickness <C 1 made 
of a composite material that periodically oscillates with period e <C 1 in in- 
plane directions. We shall derive a homogenized plate model via simultaneous 
homogenization and dimension reduction in the limit (ft,, e) — > when the ma- 
terial period e and the thickness h are either comparable or behave as e ^ ft, 
see Theorem 2.4 below. The derived model is sensitive to the relative scal- 
ing of ft and e. Our result generalizes recent results from [Neul2] where the 
one-dimensional case is studied. Regarding plates, related results have been 
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obtained for different energy scalings: In [BFFOO, BB06] the membrane regime 
has been considered. Recently, 

the energy scahng corresponding to the von-Karman plate model was studied 
in [NcuVcl]. 

This article is organized as follows. Section 2 introduces the general framework 
and discusses the main results. In Section 3 we recall the notion of two-scale con- 
vergence and characterize the two-scale limit of nonlinear strains. In Section 4 
we prove our main result. 



2 General framework and main result 

From now on, S' C denotes a bounded Lipschitz domain whose boundary is 
piecewisc C^. The piccewisc C^-condition is necessary only for the proof of the 
upper bound and can be slightly relaxed, cf. [Hollb]. 

For h > and / := (^^7 we denote by fth := S x hi the reference con- 
figuration of the thin plate of thickness h. The elastic energy per unit volume 
associated with a deformation w'' : il/j — > is given by 



y [ W{-,Vv''{z))dz. 



(3) 



Here and below z' = (zi,Z2) stands for the in-plane coordinates of a generic 
element z = (zi, Z2, 23) e Vlh and is a energy density that models the elastic 
properties of a periodic composite. 

Assumption 2.1. We assume that 

W : R2 X R3x3 ^ ^ ^r{^y^ 

is measurable and [0, l)^-periodic in y for all F. Furthermore, we assume that 
for almost every y G M^, the map M.^^^ 9 i-> W{y,F) e [0,c>d] is continuous 
and satisfies the following properties: 

(frame indifference) (FI) 
W{y, RF) = W{y, F) for ah F e R^xs^ ^ ^ SO{i); 

(non degeneracy) (ND) 
W{y,F) > cidist^(F,S0(3)) for ah F € R^^^; 

W{y,F) < C2dist2(i^,SO(3)) for all F e R^^^ ^^^^i dist2(F, S0(3)) < p; 



(quadratic expansion at identity) (QE) 
,. W{y,I + G)-Q{y,G) ^ 
G™ |G|2 

for some quadratic form Q{y, •) on R^^'^. 
Here ci , C2 and p arc positive constants which are fixed from now on. 
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We define ft := S x I. As in [FJM02] we rescale the out-of-plane coordinate: 
for X = (x'jXa) G il consider tlie scaled deformation u^{x',X'i) := v'^{x',hx3). 
Then (3) equals 

£''-'{u'') -.^ ( W{-,Vhu''[x))dx, (4) 
Jn £ 

where V ■= (V'm'', j^d^u'^) denotes the scaled gradient, and V'u'' := 
[diu'^, d2U^) denotes the gradient in the plane. 

We recall some known results on dimension reduction in the homogeneous case 
when W{y,F) — W{F). As explained in [FJM06] a hierarchy of plate models 
can be derived from := E^--^ in the zero-thickness limit /i — > 0. The different 
limiting models are distinguished by the scaling of the elastic energy relative 
to the thickness. In [LDR95] it is shown that the scaling f ^ 1 leads to 
a membrane model, which is a fully nonlinear plate model for plates without 
resistance to compression. In the regime 8^ ^ h'^ finite energy deformations 
converge to rigid deformations and, as shown in [FJM06], h~^£^ converges to a 
plate model of "von-Karman"-type. n 

In this article we study the bending regime £^ r^h'^, which, as shown in [FJM02], 
leads to Kirchlioff's nonlinear plate model: as /i — > the energy h~^£^ T- 
converges to the functional (1), with : K^^^ — > K given by the relaxation 
formula 

/ 2 

Q2{A) = minQ \ 'S^ Aapica® ep) + e^i 

\a,/3 = l 

here, Q denotes the quadratic form from (QE). 

We will see that in the non-homogeneous case the effective quadratic form Q2 
is determined by a relaxation formula that is more complicated and requires 
the solution of a corrector problem. In particular, our analysis shows that in- 
plane oscillations of the deformation couple with the behavior in the out-of-plane 
direction. As a consequence the effective behavior will depend on the relative 
scaling between the thickness h and the material period e. To make this precise 
we assume that e and h are coupled as follows: 

Assumption 2.2. Let 76 (0,oo] denote a constant which is fixed throughout 
this article. We assume that e = e{h) is a monotone function from (0,c») to 
(0, 00) such that e{h) and 7 as /i — > 0. 

The effective behavior of the homogenized plate with reduced dimension can be 
computed by means of a relaxation formula that we introduce next. We need to 
introduce some function spaces of periodic functions. From now on, Y = [0, 1)^, 
and we denote by y the set Y endowed with the torus topology, so that func- 
tions on y will be y-periodic. 

We write C{y), C''{y) and C°°{y) for the Banach spaces of y-periodic func- 
tions on that arc continuous, fc-timcs continuously differentiablc and smooth, 
respectively. Moreover, H^{I x y) denotes the closure of C°° {I , C°° (y)) with 
respect to the norm in H^{I x Y) and we write H^{I x y) for the subspace 
of functions f € H^{S xy) with Jfj^yJ^ ~ ^' '^^^ definitions extend in the 
obvious way to vector- valued functions. 
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Definition 2.3 (Relaxation formula). Let Q be as in Assumption 2.1. For 
X3 e / and A, B e M2x2^ ^^^^^ 

A{x3,A,B) := (Bap + X3Ac,i3)ie 

\a.p=i 

(a) For 7 e (0, oo) we define Q2.7 : Kyi ^ [0, 00) by 

Q2,j{A) := inf // Q(y, A{x3,A,B) + iVyq^, ^d3(b)) dydx3 

where the infimum is taken over all B € Kg-^,^ and ^ € x 3^, 

(b) For 7 = 00 we define Q2.00 ■ ~^ [0' °°) by 



Q2, 00(A) inf // Q{y, A{x3, A, B) + {Vy(t), d)) dydxs 

B.<t>4JJixY 

where the infimum is taken over all B £ Kgym, </> € i^(/, ^^^(3^, K^)) and 



KirchhofF's plate model is defined for pure bending deformations of 5 into M'^; 
precisely: 

Wl^{S,R^) := \ue W^^^{S,R^) : u satisfies (2) a.e. in S* i. (5) 



With each u € Wg''^{S) we associate its normal n := diu A d2U, and we define 
its second fundamental form 11 : 5* ^ I^sym by defining its entries as 

IIq/3 = daU ■ dpn — —dadpu ■ n. (6) 

We write II'' and n'' for the second fundamental form and normal associated 
with some G Wg''^{S, R'^). The F-limit is a functional of the form (1) trivially 
extended to L'^{n, M^) by infinity: for 7 G (0, 00] define £^ : L'^{n, R^) [0, 00], 



[ Q2,-y{ll{x'))dx' if W^f^(5,R3) 
J s 



£^{u) := { Js 

+ 00 otherwise. 

We tacitly identify functions on S with their trivial extension to = S x I: 
above u e Wg'{S,R^) means that u{x',X3) = u{x') := jju{x',z)dz for almost 
every X3 <£ /, and u e Wf'{S,R^). Our main result is the followine;: 

Theorem 2.4. Suppose that Assumptions 2.1 and 2.2 are satisfied. Then: 

(i) (Lower hound). If {m''}/i>o is a sequence with — j^^u^ dx u in 
L^{n,R^), then 

lim inf /i-^f (w'') > £^{u). 
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(a) (Upper hound). For every u G W^'^(S', R'^) there exists a sequence {u^}h>o 
with ^ u strongly in L^(r2,R^) such that 

lim h-'^£''^'''^''\u^)^EJu). 

This theorem is complemented by the following compactness result from [F JM02] , 
which in particular shows that {£'^''^'''^^}h>o is equi-coercive on L^{fl,] 



Theorem 2.5 ([FJM02, Theorem 4.1]). Suppose a sequence u'^ G H^{n,R^) 
has finite bending energy, that is 

limsup-^ / dist^(V/iM''(a;),S0(3))da; < oo. 
h^o « Jn 

Then there exists u G Wg''^ {S,M.'^) such that 

li'* - -f dx ^ u, strongly in L^(f2,M^), 

Jn 

\JhU^ ( V'lt, n) strongly m L'^{n,R^''^), 

as /t. — > after passing to subsequences and extending u and n trivially to O. 

Theorem 2.4 and Theorem 2.5 imply by standard arguments from the theory of 
F-convergence that minimizers of functionals of the form 

j^-2gKe(h)(^,-^ + "(rescaled) dead loads", 
subject to certain boundary conditions, converge to minimizers of 

E-yi') + "dead loads", 

subject to certain boundary conditions. For details see [FJM06]. 

In the special case when W{y, F) = W{F) is homogeneous Theorem 2.4 reduces 
to the result in [F JM02] . The proof of our main result emulates their argument 
as far as possible. 

We now explain our approach. The bending regime is a borderline case in the 
hierarchy of plate models. On one hand it allows for large deformations, on the 
other hand it corresponds to small strains: By Theorem 2.5 a sequence {u^}h>o 
with finite bending energy in general converges to a non-trivial deformation. 
However, the associated non-linear strain {\/ hU^Y (SI hU^) — I converges to 
zero. Indeed, let 

^ ■= h 

denote the scaled non-linear strain associated with u''. Then due to the elemen- 



tary inequality 



VWf- 



< dist(i^,S0(3)) we find that {E'^jh^o is bounded 
in when {u'^Jhyo has finite bending energy. 

The smallness of the nonlinear strain is crucial for our extension to simultaneous 
homogenization and dimension reduction: By (QE) the elastic energy is related 
to the nonlinear strain in a quadratic way - indeed, we formally have 

l^E'^^-iu'^)^ j^Q(^,E\x))dx. (8) 
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Heuristically, the right-hand side is obtained by linearizing the stress-strain 
relation, while preserving the geometric non-linearity. 

Due to the convexity of the right-hand side in (8) only oscillations of {i?''}/i>o 
that emerge precisely on scale e are relevant for homogenization. A tool to 
describe such oscillations is two-scale convergence. In Section 3 we characterize 
(partially) the possible two-scale limits of {i?''}/!.>o- This is the main ingredient 
for the lower bound in Theorem 2.4. 

Assume converges to some bending deformation with second fundamental 
form II. Then any two-scale accumulation point of {£'''}/t>o can be written in 
the form 

X3 I "^'''^ I +E{x,y) (9) 

V 

where E : il x Y ^ ^sym a relaxation field that captures oscillations and is 
a priori "unknown". In Proposition 3.2 we prove that E has to be of specific 
form. The F-limit of h"^£'^ is then obtained by relaxation: 

i|f J^J^Q [y, X3 ^ + E{x, y) ^ dy dx, 

where the infimum is taken over all E of the specific form given in Proposition 
3.2. 

We conclude this section by discussing the dependency of our limiting model 
on the parameter 7, which describes the relative scaling between h and e. The 
relaxed quadratic form (52.7 continuously depends on 7. In fact, with [NeuVel, 
Lemma 5.2] at hand one can easily identify the limits 

limQ2,-y{A) and lim Q2,-y{A) {AeM.^''^), 

7—^0 ' 7— >oo 

which yield proper quadratic forms on R^^^ that vanish on skew-symmetric 
matrices and are positive definite on symmetric matrices. In particular, the 
limit for 7—5-00 coincides with Q2.oo- The limit for 7 — > can be identified 
as well: We introduce the dimension reduced quadratic form Q2{y-,A) for all 
A e R^''^ via 

(52(2/, A) = min Q Aap{ea <S) ep) + d(^ 63 

Then Q2,7(^) converges for 7 to 

Q2.o{A) ■■= inf // Q2 {y, A + X3B + sym(Vj,C + xsVlip)) dy dxs 

where the infimum is taken over all B e Rsy„i, C G H^{y,R^) and <^ G H^{y). 

A similar behavior has been observed in [Neul2, NcuVel] where also the case 
7 = is considered (for rods and von-Karman plates, respectively). In the 
von-Karman case (see [NeuVel]) it turns out that in the regime h ^ e{h) the 
limit 7 — > of the quadratic energy density indeed recovers the energy density 
obtained via P-convergence. It is not clear whether or not this picture extends 
to the bending regime. 
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3 Two-scale limits of the nonlinear strain 



Two-scale convergence was introduced in [Ngu89, A1192] and has been exten- 
sively applied to various problems in homogcnization. In this article we work 
with the following variant of two-scale convergence which is adapted to dimen- 
sion reduction. 

Definition 3.1 (two-scale convergence). We say a bounded sequence {f^}h>o 
in two-scale converges to / e i^(ri x Y) and we write /, if 

lim / // f{x,y)i'{x,y)dydx 

h^oJn JJnxY 

for aU V € Cg°{n,C{y)). When \\f''\\L^n) -> WfWL^nxY) in addition, we say 

that strongly two-scale converges to / and write /. For vector- valued 

functions, two-scale convergence is defined componentwise. 

Since we identify functions on S with their trivial extension to fi, the definition 
above contains the standard notion of two-scale convergence on S* x y as a 

special case. Indeed, when {f'^}h>o is a sequence in i^(S'), then — ^ / is 
equivalent to 

hm / f'ix'Mx', ^) dx' = [f fix', y)i^{x', y) dy dx' 
for ah V e C^{S,C{y)). 

The main ingredient in the proof of the lower bound part of Theorem 2.4 is the 
following characterization of the possible two-scale limits of nonlinear strains. 

Proposition 3.2. Let {u'^Jhyo be a sequence of deformations with finite bending 
energy, let u G W^'^(S', R"^) with second fundamental form II, and assume that 



w'' - -f v!^ dx ^ u strongly in L'^{il, M^), 

'n 



:= ^ ^ weakly two-scale 

for some E £ L^{n x Y;M.^''^). 

(a) If-fe {0,oo) then there exist B e L^{S,Rl^i), and (f) e L^{S, H^Ixy,S:^)) 
such that 

E(x, y)=l "^'''^ + ^^"^'^ 1 + sym hy^{x, y) , i^g^lx, y)) . 

\ y 

(10) 

(h) //7 = oo then there exist B e L^{S,Rl^^), (p £ L^{n,H\y,R^)), and 
d e L2(0,R3) with 

£;(z,y)= ( + I +sym(V,0(x,y),d(x)). (11) 
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Remarks. (i) In [FJM02] a coarser characterization of possible weak limits of 
{E'^}h>o was obtained: In the situation of the previous proposition let 
{E'^Y denote the 2x2 matrix obtained from E'^ by deleting the third row 
and column. Then it was shown in [FJM02] that 

{E'^y ^ X3ll{x') + B'{x') weakly in L''{Vt,M^,^;^)). 

Proposition 3.2 refines this by capturing, in addition, oscillations on the 
scale e. 

(ii) Proposition 3.2 still only yields an incomplete characterization of the pos- 
sible structure of the two-scale limiting strain E: it is not true that every 
E in the form of (10) (resp. (11)) can be recovered as a two-scale limit of 
a sequence of nonlinear strains. For instance, when u is affine, i. e. II = 0, 
then not every two-scale limiting strain of the form (10) with B arbitrary 
and (j) = Q can emerge. 

In our construction of recovery sequences a special role is played by the 
matrix B, which is "recovered" by corrections of the isomctry of order h. 
More precisely, these corrections are obtained by solving the equation 

B = symV'g + all for 5 : 5 ^ M'' and a : S* ^ K. (12) 

As shown in [Sch07, Lemma 3.3], equation (12) can be solved locally on re- 
gions where II 7^ 0, provided that u is smooth. On the level of these "order 
h corrections" the deformed plate behaves like a shell and the condition 
II ^ corresponds to the property that the shell is developable without 
affine region. 

An important observation is that, in spite of not giving an exhaustive 
characterization of limiting strains, the result of Proposition 3.2 is just 
sharp enough to obtain the optimal lower bound for h^'^£'^'^^'^\ This 
is because on regions where II = 0, corrections associated to B can be 
ignored, since they do not reduce the energy (as Q{y,F) is minimal for 
F = 0). 

In contrast, for rods and von-Karman plates, exhaustive characterizations 
were obtained in [Neul2, Theorem 3.5] and [NeuVel, Proposition 3.3]. 

(iii) A key technical ingredient in the proof of Proposition 3.2 is Lemma 3.8 
below. It allows us to work with pieccwise constant S'0(3)-valued approx- 
imations of the deformation gradient, as opposed to smooth S'0(3)-valued 
approximations. The latter were used in the proof of the Id case given 
in [Neul2]. In the 2d case studied here, the use of such a smooth 5*0(3)- 
valued approximation would require small limiting energy, cf. [FJM06, 
Remark 5]. Thanks to Lemma 3.8, our result is not restricted to small 
limiting energy. Incidentally, the use of this lemma also simplifies the 
proof of the convergence statement in the Id case. 

The starting point of the proof of the previous Proposition is [FJM06, Theorem 
6], which we combine with the last remark in [FJM02, Section 3] in order to 
allow for 7o < 1. 
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Lemma 3.3. Let 70 G (0,1] and let h,S > with 'fo l£ j < There exists 
a constant C , depending only on S and 70, such that the following is true: if 
u G H^iQ.^'E?) then there exists a map R : S ^ 5*0(3) which is piecewise 
constant on each cube x + SY with x G (5Z^ and there exists R £ H^(S,R'^^^) 
such that 

llV^u - + \\R- RWh^s) + h'W^mhis) < C|l dist(V^w, S0(3))|li.(o)- 

Let us recall some well-known properties of two-scale convergence. We refer to 
[A1192, Vis06, MT07] for proofs in the standard two-scale setting and to [NeulO] 
for the easy adaption to the notion of two-scale convergence considered here. 

Lemma 3.4. (i) Any sequence that is bounded in L^(f2) admits a two-scale 
convergent subsequence. 

(li) Let f e xY) and let f^ G L^{Vl) be such that f^ ^ J. Then 

^ /y f{-^y)dy weakly m L'^{n). 

(ni) Let /° and f'^ G L^{Q.) he such that f^ /" weakly in L^{n). Then (after 
passing to subsequences) we have f^ f'^{x) + f for some f G L^(flxY) 
with Jy f{-,y) dy = almost everywhere in S. 

(iv) Let /" and G H^{VL) be such that f'^ f° strongly in L^{VL). Then 
f^ f^ , where we extend trivially to Vl xY . 

(v) Let /o and f^ G ^^-5") be such that f^' f° weakly in H^iS). Then 
(after passing to subsequences) 

for some (f) G L^{S,H^{y)). 

The following lemma is an immediate consequence of Lemma 3.4, cf. [NculO, 
Theorem 6.3.3] for a proof. 

Lemma 3.5. Let li" and G £r"'^(rJ,M'^) be such that weakly in 

(a) If ^ e (0,00) then there exists 4> G L'^{S,H^{I x y,R^)) such that (after 
passing to subsequences) 

V„ii''-^(VV,O) + (V,0, ^ds^). 

(b) If-f = 00 then there exist G L'^{S,H^{I x y,K.^)) and d G L'^{n,m.^) such 
that (after passing to subsequences) 

Vhu'' ^ (VV, 0) + (V,y</), d). 
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At several places in our proof we will need to make sense of a two-scale limit 
for sequences which might be unbounded in L^, but which nevertheless have 
controlled oscillations on the scale e. In order to capture these oscillations, we 
'renormalize' the sequence by throwing away the (divergent) part which does 
not oscillate on the scale e. (For bounded sequences, this latter part gives rise 
to the weak limit, but the point here is that our sequences may be unbounded.) 
Equivalently, we weaken the notion of two-scale convergence by restricting the 
admissible test functions to functions with vanishing cell average. 

More precisely, for a sequence {j^^h>a C L'^iVt) and / G L^(0 x y) with 
Jy fi'yV) dy — almost everywhere in fl, we write 

if 

h^^Jn ^ ' JJnxY 

for all G C5"(f^) and .g e C°°(3^) with J gdy = 0. (13) 

The proof of the following lemma is straightforward. 

Lemma 3.6. Let f and e L'^{n) be such that /° weakly m L'^{n) 

and J Then — ^ f'^ + f weakly two-scale. 



For the proof of Proposition 3.2 we have to identify the oscillatory part of two- 
scale limits for renormalized functions of the form j^f^ where is either 
a sequence bounded in H^{S) or piecewise afRne with respect to the lattice 
£{h)l? . The following two lemmas treat these situations. 

Lemma 3.7. Let f and f'' G H^{S) be such that f weakly in H^{S) 

and assume that 

fh 2,7, yy/ ^0 



for some 4> G L^iS, H^{y)) with fy 0('7 2/) dy = almost everywhere in S. Th 



en 



f 



h 



OSC.'J 



e{h) 

Proof. Since is independent of X3, we must show that 

-^J f\x')g{^mx')dx' ^ [[ <p{x',y)g{yMx')dx' (14) 

for all g £ C°°{y) with Jy g dy ~ and ip £ C^{S). For simplicity we write e 
instead of e{h). Let G denote the unique solution in C^{y) to 



-AyG^g, J Gdy = 
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Set Cix') := eG(^) so that 

Hence, the right-hand side of (14) equals 

= j f'(^A{G^ij)-2VG^ -ViP-G^Ai'^dx' (15) 

= - /" V/''-V(G''V)da;'-2 / f\VG^' -Vij)dx'~ [ f'^G^A^^dx'. 
J s Js Js 

By the chain rule and the definition of G'' we have 

V(GV)(-^') = VG''(x')V'(^')+G''(a;')'/'(a;') 
= V,G(|^),/'(x') + £G(^)^(x')- 
Since the right-hand side strongly two-scale converges to yG{y)ip{x), and be- 
cause V/'' V/°(x) -I- Vy4>{x,y) by assumption, we deduce that 

V/'' • V(G''V)rfa;' 

Vf{x')+Vy<P{x\y)) ■ (VyG{y)i;{x'))dydx' 

(Pix',y)AyG{yMx')dydx' 

SxY 

Hx',y)9{y)il^{x')dydx'. 

SxY 

Hence it suffices to show that the second and third integral on the right-hand 
side of (15) vanish for /i — >• 0. We treat the second integral. Since VG''(x') = 
VyG(^) strongly two-scale converges to \7yG{y), and because f'^Vtp — >■ f^S/ip 
strongly in L'^{S), wc deduce that 

-2 [ /''(VG'' • VV')da;' -> - [[ f {x')\/ij{x') -X/ yG{y) dy dx' = Q. 

Js JJ SxY 

The third integral on the right-hand side of (15) vanishes simply because f^Aip 
is bounded in ^2(5-) and G^ in L^{S). □ 

Lemma 3.8. Let f and e L°^{S) he such that ^ /° weakly-* in L°°{S). 
Assume that £ L°°{S) is constant on each cube x+e{h)Y, x G £{h)l? . Then 
we have 

^^/'^(x')V(x')ff(^) dx' ^ jj%x')Vi^{x') dx' ■ j^g[y)y dy (16) 



SxY 



for all g G C{y) with Jy 9 = and i' E G^{S). In particular, if f° S W^'^{S) 
we have 



e{hy 
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Here we write 



Proof. We first argue that (16) combined with /° g W^'°°{S) implies the con- 
vergence of j^f^ ■ Indeed, since is independent of X3 it suffices to consider 
test functions g and -0 as in identity (16). Now the statement simply follows 
from the observation that the right-hand side of (16) becomes 

{y ■ V')f{x')i;ix')g{y) dy dx' 

SxY 

by an integration by parts. We prove (16). For simplicity we write e instead of 
e{h). We denote by V''' an approximation of ip that is constant on each of the 
cubes + eY, ^ e eZ^, say 'ip^{x) := i^{£,x) where S,x & denotes the cube 
+ in which x lies. Then we have 

t^^ix)g ) dx = / /^x)^^:M^, (f ) dx (17) 



because 



dx = 0, 



since g has zero average over Y, and and ip'^ are both piecewise constant. 
Let us compute the right-hand side of (17). As £^x G eZ^ and g G C{y), we have 

and sec (after extending f^' to by zero) 

= E /'(^) / f + ^^)) • 7) (7) 



= E /"(^)^' /" f + ^-^2^) • 2^) 9{y) dy 

E /'(O^' /" vV(0-2/.9(y) rfy. 

JY 
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The first term on the right-hand side converges to zero as ft. — > because 

|VV(e + tey)-VV'(e)l <Ce 

for all t £ [0, 1], simply because V''0 is Lipschitz. 
Hence it remains to compute: 

s^f\0^'m ■ f ygiy) dy 

= E f/'(OVV(0 - / /"(^)VV'(^) dz 
+ V £2 i /^(z)VV(;^) dz ■ [ yg{y) dy 
= E f/"(C)vV(0 - / /'(^)VV(z) dz) ■ [ 



ygiy) dy 



ygiy) dy 



+ / f''ix)V'^ix) dx ■ / yg{y) dy. 

Since spt tjj G S, the last term equals 

fix)Vi;{x) dx- f ygiy) dy. 
Jy 

The claim follows because /° and because 

E f/'(C)VV(0 - / /"(^)VV(z) dz) ■ [ 

f^-^d V J£+eY I Jy 



ygiy) dy 



as ft, 0. To sec this, wc compute recalling that /''(x) = /''(^) for all x e ^+eY: 

f^^irn'm - i /'(^)vv(^) dz = f\o (vm - i vv(z) dz) < ch, 

Ji+eY \ J^+eY / 

again because VV is Lipschitz. 

□ 

Proof of Proposition 3.2, case 7 G (0, cxj). Step 1. Without loss of generality 
we assume that all u'' have average zero. Theorem 2.5 then implies that 

Vhu'' ^ R:= iV'u, n) strongly in ^^(f^^ ]r3x3) (^g-j 

where n denotes the normal to u. Let R'\ i?'' be the maps obtained by applying 
Lemma 3.3 to u'^ with (5(ft) = e{h). Due to the uniform bound on V' R^ given by 
Lemma 3.3, R^ and R^ are prccompact in L^(S', R^^'^). Hence, (18) combined 
with \\R^ — V/iu''||i2 — )• (which also follows from Lemma 3.3) shows that i?'' 
and R^ strongly converge in iF'iS, K^^^) to i?. Following [FJM02], we introduce 
the approximate strain 

, {R^YVhu\x)-I 
G [x) = . (19) 
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We set u^(a;') = JjU^^ix^xs) dxs and define z'* € H^{n,M.^) via 

u''{x', X3) = u^{x') + hx3R^{x')e3 + /iz''(a;', 2:3). (20) 
Then clearly jj z''(a;', 3:3)^x3 = and we compute 



For a given matrix A/ e M^^-^, we denote by M' the 3 x 2-matrix obtained by 
deleting the third column. We use the notation [y ■ V')R{x') ;= yidiR{x') + 
V2d2R{x') 

Step 2. Let us for the moment take for granted that there exist B' e 
L2(S',K3x2)^ ~ g L-2{^S,H^{I X 3^,R3)), V, w e i2(5^_ffi(j;^R3)) and w° e 

L^(S', R'^), such that, after passing to a subsequence, 

V/.z" ^ (VyF, ^^32), (22) 

7 

7 ^ i?'(x-') + -(yV')i?'(x-') + V,?y(x',2/), (23) 

/i 7 

a.3V'i?''e3 ^ X3V'i?(x')e3+X3Vyw(x',y), (24) 

i(i?''e3-i?''e3) ^ i(2/-V')i?(x')e3 + i^°(x') + -^J(x',y). (25) 
/i 7 7 

We now proceed to prove that the proposition follows from these convergences. 
First notice that it suffices to identify the symmetric part of the two-scale limit 
G of the sequence . Indeed, since \/ {I + hFY{I + hF) ^ I + h sym F up to 

terms of higher order, the convergence G — — ^ G implies E — symG (see e.g. 
[Neul2, Lemma 4.4] for a proof). 

We now identify symG. By combining (22) - (25) with identity (21), we find 
that i?''G'' weakly two-scale converges to 

(S', 0) + ^ X3 V'i?(x')e3 ' ^ + Vy0, ^^3^ ^ + [y ■ V')i?(x') (26) 

where 

(/)(x,2/) := z{x,y) + v{x',y) + X3w{x' ,y) + X3W° {x'). 

Due to the strong i^-convergence R^ R, we deduce that G'' weakly two-scale 
converges to (26) multiplied with i?* from the left. The first and second term 
yield 

Bix') + X3ll{x') [j 
6i(x') b2{x') 

where B denotes the 2 x 2-matrix obtained by deleting the third column of 
R^B' and (fei, 62) are defined as the entries of the third row of R^B' . Upon left 
multiplication by i?*, the last term in (26) yields a skew-symmetric term. Thus 
we have shown: 

symG(.,y) = I «y-^ + -3n I ] ^^^^^ f ^^^^ 1^^^ 







7 
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where 

bi{x') 

{x,y) := R\x')(t>{x,y) +-fX3 ( ^2(2;') 



Step 3. It remains to prove (22) - (25). Since V/jz'' is uniformly bounded in 
and since Jj z'' dxs = by construction, (22) directly follows from Lemma 

3.5. 

Next we prove (23). By Lemma 3.3 and Lemma 3.4 (i) there exists V G L'^{S x 
y, M'^^^) such that (after taking subsequences) 

Y:^_m. K (27) 

Let us verify that 

V{x', y) = B'ix') + {y ■ V')R'{x') + Vyv{x\ y), 

where B' := Jy V{-,y)dy and v £ L'^{S, {y)). This is equivalent to showing 
that 

V{x',y):{V^G){yMx')dydx' 

SxY 

{yV')R'{x'):{V^G){y)i^{x')dydx' (28) 

'Sxy 

for aU G € G^{y,m?), V € G^{S). (Here and below V^^ := {~dy.,,dy,)). Set 
^), so that {VTG'^ix') = (V^G)(j^; 



G\x') :=e{h)Gi^), so that (V')^G''(x') = (V^G)(j^). 
To prove (28), note that since / V'u'* : V'^(G''V) = 0, we have 



V'u'^ix') .y,ic'"(x')V^(x') dx' 



h 



s h 



G''{x')®\7' \p{x') dx' 



e{h) 



V'u" : G(^) ® y'^Mx') dx'. 
s ' 



The right-hand side converges to 0, since ^j^Vu'^ strongly converges in and 
G(-7^) ^ weakly in L^. In addition, Lemma 3.8 yields 



R 



and thus 



f V'u\x')-{R'^y{x') , / x' \ 



ipix') dx' 



-iyV)R'{x'):S/jGiy)^Pix') dydx'. 

SxY 7 
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On the other hand, the left-hand side converges to 

: V^G(y) V(x') dy dx' . 

I SxY 

Hence, (28) and thus (23) foUows. 

We prove (24) and (25). By Lemma 3.3 the right-hand side in (25) is uniformly 
bounded in L'^{S,M?) and thus we have (after passing to subsequences) 



[R'^ - R'')e3 vr^ 
h ' 



w{x',y) 



for some w e L'^{S x y, M^)^ get w"(.t') fyw{x',y)dy. Since R'^e^ ^ 
Rej, weakly in _ff^(S', M^), we know from Lemma 3.4 (v) that there exists w £ 
L^(S,H^{I X y,R3)) such that 

V'R''e3 ^ V'Re3 + VyW. (30) 

This implies (24). The combination of (30) with Lemma 3.7 yields ^ osc,-i^ 
^~^w. Together with (29) we get 

'^^-^w{x\y) + -i-\yV)R{x')e3, 



{R^ - R^) 63 osc,7 

h 

and (25) follows from Lemma 3.6 



□ 



Proof of Proposition 3.2, case 7 = 00. The argument is similar to the case 7 S 
(0,cxj). Therefore we only indicate the required modifications. Step 1 (of the 
proof for 7 € (0,oo)) holds verbatim modulo the following change: As a dif- 
ference to 7 G (0,03), in the case 7 = 00 we set 5{h) := [j^l^C*) where [s] 
denotes the smallest, positive integer larger or equal to s. By construction 6{h) 
is an integer multiple of e{h) and we have ^ 1. Hence, Lemma 3.3 yields 

maps R^ and i?'' with bounds uniform in h, and moreover i?'' is constant on 
each cube x + e{h)Y, x € £{h)'I? . 

Similar to Step 2 (of the proof for 7 e (0, 00)) the statement of the proposition 
can be reduced to show that (up to subsequences) 

(31) 
(32) 

(33) 
(34) 

where R is defined as in the case 7 e (0,oo), and B' € L^(S',R^^^), z € 
L'^{n,H^{y,R^)), d' e L^{n,R^),v, w e L^iS,H^{y,R^)) and w° e L^{S, 
Indeed, by the same arguments as for 7 £ (0, 00), (31) - (34) imply that 

symG(x,y)= ( '^"'^ + ''3" | + sym (V,0(x, y) , d(x)) , 






2,7, 


{VyZ^d'), 


h 


2,7^ 


B'ix') + Vyvix',y), 


xaVR'^es 


2,7, 


X3V'R{x')e3 + X3Vyw{x', y), 


^(R'^e.-R'^e,) 


2,7, 


w"{x'). 
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where 

(j){x,y) -.^^ R\I+v + X3w), d^R*d' + R*w° 



and B, R and {bi, 62) are defined as in tlie case 7 € (0, 00). 

The proof of (31) - (34) is similar to Step 3 of the proof for the case 7 € 
(0,00). □ 




4 Proof of Theorem 2.4 

As a preliminary step we need to establish some continuity properties of the 
quadratic form appearing in (QE) and its relaxed version introduced in Defini- 
tion 2.3. 

Lemma 4.1. Let W be as in Assumption 2.1 and let Q be the quadratic form 
associated to W through the expansion (QE). Then 

(Ql) Q{-,G) is Y -periodic and measurable for all G € R^^'^, 

(Q2) for almost every y G the map Q{y, •) is quadratic and satisfies 

cil symG|2 < Q{y,G) = Q{y,symG) < C2I symCj^ for all G e R'''^^ 

Furthermore, there exists a monotone function r : M+ — > U {+00} that can 
be chosen only depending on the parameters ci,C2 and p, such that r{S) — > as 
(5—^0 and 

VG G R3^3 . / + G) - Q{y, G)| < |Gpr(|G|) (35) 

for almost every ?/ G M^. 
(For a proof see [Ncul2, Lemma 2.7].) 

Lemma 4.2. (a) Let 7 G (0,oo). For all A G Rgyi^ there exist a unique pair 
{B, (j)) with B G and (pG H^L xy, M^) such that: 

Q2,j{A) ^ jj^ ^Q(y, A(x3, A, B) + (V,0 , ^d^^) ) dydx^ 

The induced mapping R^y^ 3 A H' (B, 0) G Mgy^ xiJ^(/xy,]R^) is bounded 
and linear. 

(b) Let 7 ~ 00. For all A G ^sym there exist a unique triple {B,4>,d) with 
B G R"^^^,, (f) G L'^{I,H^{y,m.^)) and d G L2(/,R3) such that: 

Q2.oo(A) ^11 Qiy, A{x3, A, B) + (Vj,<^ , d) ) dydxs 

JJixY 

The induced mapping ]R2x2 g ^ ^ {B,(f>,d) G M2x2 ^ L'^{I, H^{y,M.^)) x 
X^(/, R^) is bounded and linear. 
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Proof of Lemma ^.2. The proof is standard. 

We only comment on the proof of part (a), i. e. for 7 G (0, 00). Wc start with 
the following Korn-type inequality 



V^ei/i(/x3;,M3) : |lVil,.(,,y,R3) <Cyy^^^|sym(V,^ , ^d^^ (36) 

for some constant C > 0. Since Q is elliptic in the sense of Lemma 4.1 (Q2), 
and since Q vanishes for skew-symmetric matrices, for each pair (A, B) of sym- 
metric 2x2 matrices we can find a unique function = <\>a.b G ll^ij x y , R-^) 
minimizing the integral 



Q^{A,B) (y, A(x3,A,B) + (V^(/. , i^s'/') ) dydx^. (37) 



Evidently (/) depends linearly on (A, i?). In particular, for i = 1,...,3 there 
exist E, GH^ilx y,R^^i), F, e H\l x y,R^^^) such that 

Mx3,y) = E,{x3,y) : A + F,{x3,y) : B. (38) 

As a consequence Q-y is a quadratic form and it is easy to check that there exist 
positive constants c-y^i < c-y^2 such that 

c^,i(|^P + \B\') < Q^{A,B) < c^,2(|Ap + \B\'). 



Hence, we conclude that there exists a bounded, positive definite operator 

X R2X2 g^^h that 

Q^iA,B) = {A{A,B), (symA,B)), 



A : R2X^ X R2X2 ]R2X2 ^ r2x2 ^^^^ that 



where (•, •) denotes the inner product in Ml^y^^ x Mg^J^. By a block- wise decom- 
position of the right-hand side wc get 

g(A, B) = AiA : A + A2{A) : B + A3B : B 

where Ai,A2,A3 E L(Mgyjjj) arc bounded operators on ]Rsy„j and Ai and ^3 
arc positive definite. A straightforward calculation shows that this expression is 
minimized (with respect to B) by Ba — ~^A3^A2 sym^. Since this expression 
is obviously linear in A, the desired pair of functions associated with A is given 
by Ba and </) = (f>A,BA ■ ^ 



Proof of Theorem 2.4 (lower bound). Without loss of generality we may assume 
that j-^u^dx = and limsup^_jg ''■^^''-'(u'') < 00. Wc only consider the 
case 7 G (0,oo). The proof in the case j = 00 holds verbatim. In view of 
(ND), the sequence has finite bending energy and the sequence E'\ see 
(7), is bounded in L^(fi, R'^^^). Hence, from Theorem 2.5 we deduce that 
u € Wg''^(S,M.^). By Lemma 3.4 (i) and Proposition 3.2 (i) we can pass to 
a subsequence such that, for some E G i^(fi x y;R^^^), 

E'^ ^ E, 
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where E can be written in the form of (10). As in [FJM02] a careful Taylor 
expansion of W{^^,I + hE^{x)) combined with the lower semi-continuity of 
convex integral functionals with respect to weak two-scale convergence (see e.g. 
[Vis07, Proposition 1.3]) yields the lower bound 

\immi ^£^^<^\u'') > [[ Q{x\y,E{x,y))dydx. 
Combined with (10) the right-hand side becomes 




where we used that Q{x' ,y, F) only depends on F. Minimization over B G 
L2(S',M2x2) and 0e L^{S,H\I x y,R^)) yields 

liminf-^£''^"('')(w'') > / Q2,^{ll(x'))dx' ^S^iu). 

□ 

It remains to prove the upper bound. As in [Sch07] and other related results, the 
key ingredient here is the density result for W'^'^ isometric immersions estab- 
lished in [Holla, Hoi lb] (cf. also [Pa04] for an earlier result in this direction). 
It is the need for the results in [Hoi lb] that forces us to consider domains S 
which are not only Lipschitz but also piecewise - more precisely, we only 
need that the outer unit normal be continuous away from a subset of dS with 
length zero. 

For brevity, we denote by AiS) the set of all u e W^'^iS, R^) n C°°(S, R^) with 
the following property: 

For each B e C°"(S,R'^^^) with satisfying B = in a neighborhood of {x G 5 : 
U{x) = 0}, there exist a e C°°(S) and g e C°°(S,R^) such that 

S = symV'.g + aH. (39) 

The key ingredient in the proof of the upper bound is the following lemma. 

Lemma 4.3. The set A{S) is dense in W^'^(S') with respect to the strong W^'^- 
topology. 

Proof. This follows by combining the construction from [Holla, Hollb] with 
the arguments in [Sch07] leading to his Lemma 3.3. His result was recently 
re-derived in [HLP] in a slightly different context. □ 

Thanks to Lemma 4.3 it will be enough to construct recovery sequences for lim- 
iting deformations belonging to A{S). First we present a construction assuming 
the existence of a and g satisfying (39). 

Lemma 4.4. For u G Wf^{S,R^) n W^-°°{S,R^), a e W^^°°{S) and g e 
M^1'°°(5',M2) define 

B{x') = sym V'g(x') + a{x') ll{x'). 
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(a) Letj e {0,00) and(j)€ C^iS,C°°{Ixy,R^)). Then there exists a sequence 
W'}h>a C i/i(rj,R3) such that 



V ^ (Vm, n) uniformly in Q. 



lim ^£^''^^\u'') 



Q (y, A(II(a;'),B(x')) + (Vy^, ^83^) ) dydx^dx'. (40) 



(b) Let = 00, (j) & Cf{n,C°°{y,B?)) and d G C^{Vt,M?). Then there exists 
a sequence {u^}h>o C iJ^(51,R'^) such that 



(V'm, n) uniformly in £7 



lim 



Qiy, AiUix'),B{x')) + iVycj), d)) dydx^dx'. (41) 



Proof. We start with the case 7 S (0, 00) and foUow [Sch07, Theorem 3.2 (ii)]. 
Consider 

i;''(x-) := u{x') + h\{x:i + a{x'))n{x') + {g{x')-V')u{x') 
R{x') (V'u(.t'), n{x')). 

Here (5- V')u stands for J2a=i 2 9o'^aU. A direct calculation (sec [Sch07]) shows 
that 

R*Vhv'' =I + hG 



with 



G 



X3 II(x') + B 



B:=V'g + aU, b := - 



dia 
d2a 



llg. 



To v'^ we add an oscillating correction and a term compensating for 6: 
u'^ix) := z;''(x)+te(/i)0(x,^) 



(42) 



where 



(x,?/) i?(x-') (^(x, ?/) + 7x3 



Since a, g and (/> are sufficiently smooth, the uniform convergence of and 
V hU^ directly follows from the construction. Moreover, we have 



i?*Vhw'' = / + h{G + G^) 
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where 



Here the remainder term C satisfies hmsup^_^o < oo. Again, since u, <j>, 
a and g are sufficiently smooth, we have 



limsup/isup (|G(x)| + |G''(x)|) = 0. 



(43) 



Hence, (FI), (QE) and (35) yield 



lim sup 



}_^H.e(H) _ Q(_^^ G(X) + G\x)) dx 



0, 



and it suffices to show that 



1™ I Q{^,G{x)+G^{x))dx=[YL. H. S. of (40)]. 



(44) 



By construction the sequence G + strongly two-scale converges to some limit 



symG = A(H,symB) + sym ( \/y(t){x,y), \dz4){x,y) 



Hence, by the continuity of convex integral functionals with respect to strong 
two-scale convergence we can pass to the limit in (44) (e. g. see [Neul2, 
Lemma 4.8]). Since Q{y,F) only depends on the symmetric part of F, this 
completes the proof for the case 7 € (0, 00). 

The proof for 7 = 00 is similar to the above reasoning. Essentially, we only need 
to replace (42) by 



u''(x) = v^{x) + he{h)(f>{x, j2^) + h^d{x), 

where 

{x,y) = R{x')4>{x,y) and d{x) — R{x') j d{x',t)dt + X3 

-1/2 



b{x') 




□ 



Proof of Theorem 2.4 (Upper bound). We only consider the case 7 € (0,oo) 

since the argument for 7 = 00 is the same. We may assume that £j{u) < cxj, so 

u G VF|'^(S', R"^). Moreover, since Q2.-y (see Lemma 4.2) is continuous, it suffices 

to prove the statement for a dense subclass of Wj'^(S', M.^). Hence, by virtue of 

Lemma 4.3, we may assume without loss of generality that u € ,4(5'). 

By Lemma 4.2 there exist B € L2(S',R2x2) j^^^j ^ ^ L^{S,H^iI x y,R^)) such 

that 

£^i^)= ff Q{y,A(Ll,B) + sym{Vy(P, ^d3cf>))dydx. (45) 
JJnxY 
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Since B{x') linearly depends on II(x') we know in addition that B{x') = for 
cc' G {II = 0}. 

By a density argument it suffices to show the following: There exists a doubly 
indexed sequence u^'^ € H-^{V,,E.^) such that 

linisuplimsup llu''''^ - m||i,2(j2^r3) =0, (46) 



lini sup 

/i-i-O 



< S. (47) 



Indeed, if this is the case then we obtain the desired sequence by extracting a 
diagonal sequence (e. g. by appealing to [Att84, Corollary 1.16]). 

We construct u^'^ as follows: By a density argument we may choose for each 
S>0 functions B^ G C°°(S,M.^^^) and 0* G C;?° (£•, C°° (/ x y,R^)) such that 

\\B' - B\\l2i^s) + l|Vy0* - V,0|U.(f2xy) + \\d3^' - ds^WmnxY) < 5\ (48) 
B'' = in a neighborhood of { II = }. (49) 

Because u G ^(S*, R'^) and due to (49) we can find for each 5 > smooth 
functions and such that B^ = sym V'g'' + II. We apply Lemma 4.4 to 
u, g^,a^ and (j)^ . We obtain a sequence {u'^'^}h>o that uniformly converges to 
u as ft- — > 0. Hence, (46) is satisfied. Moreover, Lemma 4.4 yields 

lim ^S^iu'^'') = ff g(y,A(II,B*) + sym(V,/, ^d^cj)') dydx. 



By continuity of the functional on the right-hand side, combined with (48) and 
(45), the bound (47) follows. □ 
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